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ABSTRACT 

« 

The  combustion  front  in  a  column  of  a  granular  mixture  of  potassium 
perchlorate  and  aluminum  has  been  observed  to  propagate  stably  at 
either  a  low  speed  (300  m/s)  with  a  short  reaction  zone  (2  cm),  or  a 
high  one  (900  m/s)  with  an  elongated  reaction  zone  (8  cm),  and  with  an 
occasional  rapid  transition  between  them.  Calculations  are  made  to 
categorize  the  observed  phenomena  as  those  of  either  a  deflagration  or 
a  detonation. 

Calculating  the  penetration  depth  from  the  front  into  the  mixture  for 
diffusion  of  molecules,  thermal  energy,  and  radiant  energy  shows  that 
these  deflagrative  mechanisms  have  negligible  effect  at  the  rates  observed 
and  suggests  that  deflagration  would  proceed  at  less  than  0. 3  m/s.  A 
calculation  of  the  deflagration  branch  of  the  Hugoniot  by  a  parametric 
technique  also  shows  that  only  rates  of  the  order  of  0. 1  m/s  exist  for 
deflagrations.  The  calculated  sound  speed  for  the  unreacted  granular 
mixture  is  less  than  300  m/s.  The  hypotheses  of  a  simple  deflagration 
and  of  a  deflagration  following  a  precompression  shock  are  rejected. 

The  RUBY  computer  program  calculated  a  detonation  velocity  for  the 
mixture  of  4600  m/s.  This  would  be  reduced  by  the  observed  lateral  rare- 
1  faction  in  the  reaction  zone.  It  is  concluded  that  to  establish  the  details 

of  the  detonatlve  nature  of  the  phenomena  further  experiments  are  re¬ 
quired  on  effect  of  confinement,  interstitial  gases,  diameter  and  length 
of  mixture  column,  and  time-resolved  spectroscopy  of  the  reaction  zone 
structure.  Future  theoretical  consideration  should  either  combine  the 
physical  kinetics  of  aluminum  particle  vaporization  with  a  two-phase 
theory  or  examine  coupling  between  shook  and  combustion  within  the  re¬ 
action  zone. 
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INTRODUCTION 


This  report  explores  the  nature  of  the  combustion  wave  in  a  granular 
mixture  of  potassium  perchlorate  and  aluminum.  The  mixture  is  a  pow> 
der,  which,  when  loose-loaded  into  a  case  and  initiated  by  an  imbedded 
small  explosive,  produces  a  short,  intense  flash  of  light.  Several  years 
ago,  research  was  started  on  the  physics  of  these  photoflash  items.  The 
first  phase  (Ref  1)  of  the  research  showed  by  means  of  flash  radiographs 
that  a  small  PETN  charge  centered  in  nonreacting  metal  dust  (simulating 
mechanical  properties  of  the  mixture)  formed  an  expanding  bubble  of 
explosion  products  surrounded  by  an  ever-thickening  spherical  shell  of 
compressed  metal  dust.  This  behavior  was  also  described  theoretically. 

The  nonreacting  metal  dust  was  then  replaced  by  the  potassium  per¬ 
chlorate/aluminum  mixture  and  its  interaction  with  the  confining  case 
was  studied  in  relation  to  the  light  output  and  the  size  of  the  cloud  that  is 
formed  from  the  item  (Ref  2) . 

Theoretical  treatment  requires  information  on  the  combustion  process 
that  occurs  in  the  confined  mixture  under  explosive  initiation.  The  next 
phase  (Refs  3,  4)  was  therefore  devoted  to  gathering  information  on  com¬ 
bustion  rates.  The  experimental  approach  used  was  to  confine  the  mix¬ 
ture  along  the  axis  of  a  thick  cylinder  of  transparent  plastic  and  then  use 
high  speed  photographic  techniques  to  observe  the  propagation  of  com¬ 
bustion  after  initiation  by  an  explosive  located  at  one  end  of  the  mixture 
column.  The  results  obtained  were  more  complicated  than  anticipated. 

In  Figure  l,  successive  sU^es  in  one  type  of  combustion  are  shown, 
as  viewed  initially  by  backlighting  and  later  only  by  emitted  light.  The 
upper  left  photo  shows  a  transparent  plastic  rod  3  inches  in  diameter  and 
20  inches  long,  with  a  3/ 8-inch-diameter  colum  of  loose-loaded  mixture 
along  the  axis.  The  notches  in  the  supporting  bar  are  1  inch  apart.  Ini¬ 
tiation  occurred  at  the  left.  The  photos  shown  are  96  microseconds  apart. 
The  lower  right  photo  is  the  stable  reaction  zone  profile  that  is  observed 
to  move  at  about  900  meters  per  second.  Note  that  the  highly  luminous 
region  extends  forward  into  the  mixture  column  and  that  there  follows  a 
dark  area  preceding  a  luminous  cloud. 

In  some  cases  a  different  stable  profile  is  observed  propagating  at 
only  about  300  meters  per  second,  with  no  such  forward  extension,  or 
dark  zone  and  following  cloud.  This  slower  profile,  a  small,  bright 


region  surrounded  by  a  much  larger,  less  luminous  cloud,  is  shown  in 
Figure  2  together  with  the  type  described  above.  In  addition  a  rapid 
transition  (almost  always  to  the  higher  speed)  between  the  two  stable 
states  was  occasionally  observed.  By  both  necessity  and  interest,  the 
research  since  then  has  been  concerned  with  e^laining  the  essential 
features  of  these  experimental  results. 

The  following  possibilities  were  previously  considered  qualitatively 
(Ref  4)  high-speed  deflagration;  deflagration  with  a  precompression  shock; 
deflagration  preceded  by  an  interstitial  detonation  (the  latter  Involves 
only  the  fines  suspended  in  the  interstices  of  the  powder) ;  or  a  detonation 
with  an  expanded  reaction  zone  associated  with  the  chemical  and  physical 
kinetics  of  the  medium.  This  report  begins  a  quantitative  treatment  by 
providing  calculations  on  diffusion  of  molecular  species  and  energy, 
sound  speeds  possible  in  the  medium,  the  deflagration  branch  of  the 
Hugoniot  curve,  and  the  Chapman-Jouguet  values  for  a  detonation.  In 
each  case  the  significance  of  the  results  is  discussed  in  determining  the 
mechanism  involved  in  the  propagation  of  the  combustion.  Directions 
for  future  e}q)erimental  and  theoretical  studies  are  then  proposed. 

It  will  be  assumed  in  the  treatment  that  follows  that  the  reader  has 
some  familiarity  with  deflagrative  and  detonative  combustion  theories 
(Ref  5,  6)  and  other  related  subjects.  However,  adequate  information 
will  be  supplied  on  the  physical  significance  of  the  equations  used,  to 
permit  following  the  logic  of  the  presentation. 


DIFFUSIVITY  CONSIDERATIONS 


A  deflagration  is  a  form  of  combustion  in  which  the  flame  progresses 
by  diffusion  of  both  activating  species  and  energy  from  the  reaction  zone 
into  the  unreacted  mixture.  The  classical  partial  differential  equation 
(Ref  7)  that  governs  this  process  is  (for  one  dimension,  constant  Jr,  no 
sources) 

77”  ~  7F  (1) 


Where  T  is  the  temperature  for  energy  flow  or  the  concentration  for  flow 
of  molecular  species.  J(  is  a  constant  called  the  diffusivity,  \diich  from 
the  equation  has  the  dimensions  of  length  squared  divided  by  time,  or  of 
velocity  times  length. 
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For  our  case,  we  replace  the  combustion  wave  which  proceeds  at  a 
constant  velocity  U,  maintaining  the  same  profile,  by  a  model  having 
energy  release  in  a  fixed  plane  at  x  =  0,  preheating  the  region  in  front 
of  this  plane  which  moves  with  speed  -U  (to  the  left)  toward  it.  At  and 
behind  the  plane  (x  =  0)  a  fixed  temperature  or  concentration  (Tq)  is 
assumed.  Although  the  model  can  be  made  more  sophisticated,  it  is  not 
necessary  to  do  so  for  the  order  of  magnitude  results  required.  We 
rewrite  equation  1  as 


(la) 


Solving  equation  la  subject  to  the  boundary  conditions  that  xi 0  T  =  T©; 
X— T-»««one  finds  for  x>0  (the  region  that  is  preheated) 


S-T- 


-'U 


and  energy  flvix  = 


(2) 


The  temperature  and  flux  decay  to  l/e  of  the  value  at  x  =  0  as  one  moves 
a  penetration  distance,  L,  to  the  right  (x  >0),  where 

U  (3) 

This  result  corresponds  to  the  dimensional  result  following  equation  1. 

The  penetration  distance  is  the  length  of  unreacted  mixture  directly  in 
front  of  the  flame  front  into  which  energy  and  molecular  species  can 
penetrate  to  initiate  chemical  reaction.  It  is  shown  later  that  the  values 
of  L  from  equation  3  are  about  1  micron. 

For  the  flame  front  to  continue  to  move  forward  with  the  same  profile 
at  a  speed  of  U,  the  energy  release  from  the  flame  front  into  the  penetra¬ 
tion  length  would  have  to  occur  in  a  time  t  =  L/U.  For  U  =  300  m/s, 

L  =  1  microm,  t  =  3.3  x  10“^  seconds  results.  However,  a  separate  cal¬ 
culation  shows  that  the  15-micron-diameter  aluminum  granule  requires 
of  the  order  of  10~7  second  to  liquefy  and  10“5  second  to  fully  vaporize 
under  optimum  conditions  of  continuous  closure  to  the  flame  temperature. 
Since  the  energy  can  not  be  released  except  in  proportion  to  aluminum 
vaporized,  one  must  conclude,  since  10“9,  that  a  deflagration 

propagated  by  these  mechanisms  at  300  m/s  cannot  exist  in  this  mixture. 


4 


tiote  that  if  particle  sizes  were  reduced  to  sub-micron  dimensions  so  that 
vaporization  were  orders  of  magnitude  faster,  the  arguments  advanced 
would  not  be  valid.  This  will  be  used  in  the  later  discussion  of  inter¬ 
stitial  effects. 

To  contrast  with  the  fact  that  a  deflagration  at  300  m/s  in  this  mixture 
is  not  possible,  consider  the  possibility  of  one  at  0.3  m/s.  Using  equa¬ 
tion  3,  it  follows  that  now  the  penetration  length,  call  it  Li  will  be  10^ 
greater  than  the  penetration  length  for  300  m/s;  thus  is  of  the  order 
of  103  micron.  Then  tj^  =  Li=  3.3  x  10~3  seconds.  These  values  of 

Ui 

and  t]^  Indicate,  since  10~3«  10~3,  that  a  reasonable  speed  for  deflagra¬ 
tion  if  it  existed  would  be  about  0.3  m/s  or  less. 

The  above  presentation  has  used  the  fact  that  L  computed  from  equa¬ 
tion  3  is  of  the  order  of  1  micron.  The  values  of  L  are  derived  below 
and  substantiate  the  conclusions  of  the  previous  paragraphs. 

Molecular  Diffusion 

Chemical  reaction  can  be  initiated  by  molecular  species,  present  in  the 
flame  zone,  moving  into  the  unreacted  mixture.  The  concentration  of 
these  activating  species  will  be  established  by  diffusion.  Typical  values 
(Ref  8)  of  diffusivity  (cmVsec)  for  gases  into  air  at  1  atm  and  about  20P  C 
are,  for  CO2,  O2,  and  H2O  vapor,  0. 139,  0. 178,  and  0. 239  respectively. 
For  dilute  gas  binary  mixtures  (Ref  9)  diffusivity  decreases  inversely 
with  pressure  and  increases  roughly  as  the  1.  8  power  of  temperature. 

The  flame  temperature  is  estimated  as  3000°  K  and  the  pressure  therein 
should  be  less  than  atmospheric.  However,  the  unreacted  mixture  is  at 
room  temperature  and  atmospheric  pressure  (except  for  the  effect  of 
simultaneous  energy  diffusion  considered  below) . 

For  estimating  purposes,  one  can  use  X=  0,2,  subject  to  a  maximum 
increase  due  to  temperature  of  63  times,  that  is,  (3000/300)^ •  ®  .  The 
lowest  observed  flame  velocity  was  3  x  10^  cm/s.  The  penetration  dis¬ 
tance  L  =  0.2/(3  X  10^)  =  0.066  micron,  which,  even  if  multiplied  by  63, 
is  only  about  4  microns. 

Thermal  Diffusion 

The  {q)proach  is  similar  to  that  above.  For  aluminum  particles  the 
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thermal  diffusivity  (Ref  8)  is  about  1.35  and  for  air  it  is  0,16.  Using  the 
larger  value  and  equation  3  we  find  L  =  0.45  micron. 

Radiant  Energy  Diffusion 


The  key  assumption  made  in  this  case  is  that  thermal  equilibrium 
exists  between  radiation  and  the  granular  media.  One  can  then  use  (Ref 
10)  the  equilibrium  (black  body)  relationships  to  show  that  the  radiative 
flux  is  related  to  the  temperature  distribution  by 


(4) 


Here  is  radiation  conductivity  analogous  to  thermal  conductivity, 
fining  a  radiation  diffusivity  s  using  equation 

find  a  penetration  distance 

L  T.  1  ^ 


so 


for  U  =  300  m/s  T  =  SSOO^K 


De- 
3  we 

(5a) 


s  for  U  =  0,  3  m/s  T  =  3500Pk  (5^) 

rC 

where  a  is  the  Stefan  Boltzmann  constant  (7,67  x  10"15  ergs/cm®  deg^), 
c  the  velocity  of  light  (3  x  lOl®  cm/sec),  ^the  density  of  air  (0. 00i25 
gms/cm^),  Cp  the  specific  heat  of  air  (0,285  cal/g),  and  the 
Rosseland  mean  free  path,  which  is  essentially  the  photon  mean  free 
path. 

Microphotographs  (Ref  1)  of  granular  media  at  this  loading  density 
and  particle  size  indicate  that  must  be  of  the  order  of  30  microns, 
that  is,  twice  the  average  particle  size.  One  can  interpret  equation  5a 
as  leading  to  a  penetration  distance  for  radiation  by  diffusion  of  only  1 
micron.  However,  it  is  more  instructive  to  compare  equation  5a  with 
5b.  For  the  latter,  the  result  states  that  radiation  will  be  absorbed  and 
re -emitted  over  a  length  30  times  the  mean  free  path  of  a  photon  {\») 
when  energy  is  transmitted  forward  by  radiation  and  both  the  flux  and  the 
local  temperature  associated  therewith  will  fall  to  1/e  of  the  maximum  of  a 
distance  L.  A  value  of  L  less  than  indicates  that  the  basic  process 
of  absorption  and  re -emission  has  no  significant  role.  This  confirms  the 
previous  numerical  analysis  for  all  modes  of  diffusion  based  on  the  1- 
micron  length. 
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It  follows  that  at  the  observed  rates  of  300  m/s  or  higher  and  maxi 
mum  temperature  of  3, 500“  K,  radiation  is  not  a  significant  factor. 


SOUND  SPEED  CONSIDERATIONS 

The  combustion  front  of  a  one-dimensional  deflagration  progresses  at 
a  rate  that  is  subsonic  with  respect  to  the  sound  velocity  in  the  unreacted 
medium  (Ref  5)  preceding  it.  For  a  deflagration  preceded  by  a  precom¬ 
pression  shock,  the  shock  is  supersonic  with  respect  to  the  medium  it  is 
entering,  and  the  deflagration  follows  at  a  subsonic  speed  measured  rela¬ 
tive  to  the  precompressed  medium  preceding  it.  A  detonation  front  pro¬ 
gresses  at  a  supersonic  rate  with  respect  to  the  unreacted  medium.  Thus 
a  comparison  of  the  sound  speed  in  the  unreacted  medium  and  that  in  a 
compressed  medium  with  the  experimentally  observed  propagation  rate 
would  establish  the  nature  of  the  phenomenon  under  the  assumption  that  it 
is  one  dimensional. 

To  consider  these  possibilities  one  requires  data  on  the  sound  speeds 
in  the  granular  medium  as  a  function  of  density.  Since  applicable  experi¬ 
mental  data  or  theoretical  results  were  not  found  in  the  literature  (Refs 
18-28)  and  experimental  facilities  for  the  measurements  were  not  avail¬ 
able,  theoretical  calculations  as  described  below  were  used.  The  num¬ 
erical  results  are  shown  in  Figure  3.  Their  significance  with  respect  to 
the  observed  combustion  phenomena  is  treated  in  the  later  subsection 
entitled  "Discussion  of  Sound  Speed  Results. " 

Assumptions  and  Definitions 


In  general,  for  a  medium  with  no  changes  in  chemical  composition, 
the  sound  speed  is  defined  by 

2 


■  &): 


(6) 


However,  different  sound  speeds  can  exist  depending  on  how  the  constant 
entropy  condition  is  met.  For  a  two-phase  (gas  and  granules)  flow  this 
will  depend  on  the  relative  motion  of  the  two  phases  and  the  heat  transfer 
between  them  during  the  period  of  the  sound  wave. 


The  equilibrium  sound  speed  is  defined  as  that  in  which  only  the  total 
entropy  is  held  constant.  Since  relative  flow  of  gas  with  respect  to 
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granules  is  a  viscous  flow,  generating  entropy,  one  must  assume  here 
that  there  is  no  significant  relative  motion  during  the  period  of  the  sound 
wave.  This  would  be  valid  for  the  frequency  low  enough  so  that  the  time 
required  to  equalize  velocities  is  a  small  fraction  of  the  period.  Further, 
since  only  the  total  entropy  is  held  constant,  heat  transfer  between  the 
two  phases  is  permissible  if  there  is  no  generation  of  entropy.  This  will 
be  the  case  if  the  transfer  is  quasistatic  (i.  e.,  takes  place  with  infinitesimal 
temperature  differences),  which  can  be  met  only  for  sound  waves  with 
periods  long  compared  to  the  time  required  for  temperature  equalization. 

To  make  quantitative  estimates  on  how  low  the  frequency  of  a  sound  wave 
must  be  for  an  equilibrium  sound  speed  to  exist,  one  would  have  to  cal¬ 
culate  the  viscous  relaxation  time  and  the  thermal  relaxation  time  for 
the  medium  and  make  the  period  several  times  greater  than  either.  For 
our  immediate  purpose,  it  is  adequate  to  think  of  the  equilibrium  sound 
speed  as  that  valid  in  the  low  frequency  limit. 

Next,  we  define  two  different  frozen  sound  speeds  in  each  of  which  we 
specify  how  the  total  entropy  is  maintained  constant  by  imposing  additional 
conditions.  Consider  that  frozen  sound  speed  (designated  I)  in  which  there 
is  no  heat  flow  and  no  relative  motion  between  the  two  phases.  The  entropy 
of  each  phase  is  maintained  constant  separately  and  no  entropy  generation 
is  permitted,  so  obviously  the  total  will  be  constant.  Physically,  no  heat 
transfer  requires  a  high  frequency  sound  wave  (short  period)  so  that  no 
significant  heat  flow  can  occur.  Velocity  equalization  has  the  opposite 
requirement  of  long  periods.  However,  one  can  create  some  special 
cases  where  both  conditions  might  be  simultaneously  satisfied.  Suppose 
the  medium  is  sufficiently  compressed  so  that  the  aluminum  granules 
form  a  lattice  with  pockets  of  gas  almost  totally  enclosed  by  metal  sur¬ 
faces.  A  rapid  compression  of  the  lattice  will  transport  the  gas  with 
very  little  relative  motion  and  for  such  short  times  the  heat  transfer 
could  be  negligible  also.  The  frozen  sound  speed  (I)  is  therefore  the  high 
frequency  limit  for  the  porous  lattice.  However,  for  the  granular  aggre¬ 
gate  model  (to  be  defined  later)  there  is  no  frozen  sound  speed  (I)  because 
there  is  always  a  generation  of  entropy  due  to  friction  as  the  granules  move 
relative  to  each  other,  so  that  one  cannot  hold  the  total  entropy  and  also 
that  of  the  gas  constant. 

Another  frozen  sound  speed  (designated  II)  would  also  permit  no  tem¬ 
perature  equalization,  but  would  have  the  particles  stationary  during  the 
gas  motion.  For  sufficiently  high  frequencies  the  wavelength  of  the  sound 
wave  is  a  fraction  of  the  particle  size  and  the  relative  motion  of  the  gas 


with  respect  to  the  particles  produces  an  insignificant  quantity  of  entropy. 
This  frozen  sound  speed  (II)  is  the  high  frequency  limit  for  the  dusty  gas 
model  and  the  granular  aggregate  model  (particles  stationary) ,  and  will 
correspond  with  the  ordinary  sound  speed  in  the  gas  alone. 

It  is  to  be  noted  that  physically  there  are  cases  where  the  total  entropy 
is  not  constant  during  the  passage  of  a  weak  mechanical  vibration  or  sound 
wave.  The  theoretical  treatment  must  therefore  calculate  both  the  equilib¬ 
rium  and  frozen  sound  speeds,  which  are  the  low  and  high  frequency  limits, 
respectively,  and  assume  that  the  intermediate  nondefinable  cases  repre¬ 
sent  velocities  between  them.  This  range  is  then  used  for  comparison  with 
the  observed  velocity  of  propagation  of  combustion  to  determine  whether 
the  latter  is  subsonic  or  supersonic. 

Next,  we  consider  assumptions  on  the  nature  of  the  medium.  In  the 
limit  of  low  densities,  the  powder  is  a  dusty  gas.  In  fact  before  this 
lower  limit  is  reached  there  may  be  small  pockets  of  powder  which  act  as 
a  dusty  gas.  This  is  because  a  dusty  gas  is  defined  as  one  in  which  the 
particles  do  not  participate  directly  in  resisting  the  ^plied  load.  Par¬ 
ticles  may  be  in  contact  but  if  they  are  not  under  load,  the  medium  is  still 
called  a  dusty  gas.  In  the  limit  of  high  densities,  the  powder  is  a  porous 
elastic  solid  with  a  lattice  formed  by  the  particles  and  an  interstitial  com¬ 
pressible  gas.  In  the  range  of  densities  between,  there  is  a  relative 
motion  of  particles  forming  successive  load-sustaining  geometrical  con¬ 
figurations.  This  intermediate  region  is  called  the  granular  aggregate 
case  and  is  believed  to  correspond  most  closely  to  the  actual  medium. 

It  will  be  treated  by  using  an  experimentally  determined  stress-strain 
curve.  The  dusty  gas  and  porous  elastic  solid  limiting  cases  will  also 
be  considered. 

In  all  cases,  the  medium  is  taken  as  consisting  of  aluminum  granules 
only  instead  of  as  a  mixture  of  aluminum  and  potassium  perchlorate.  The 
experimental  curve  used  for  the  granular  aggregate  could  only  be  safely 
determined  for  aluminum  alone.  The  equations  which  will  be  used  are 
derived  only  for  a  metal  powder,  and  extending  these  to  more  than  one 
solid  component  would  involve  considerable  algebra.  However,  this 
assumption  introduces  no  significant  error  because  the  difference  between 
the  mechanical  properties  of  air  as  compared  to  those  of  either  aluminum 
or  potassium  perchlorate  is  orders  of  magnitude  greater  than  the  difference 
between  the  mechanical  properties  of  the  two  solids.  For  this  reason  the 
conclusions  based  on  an  aluminum  dust  are  the  same  as  would  follow  from 
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an  extended  consideration.  This  will  be  evident  in  the  numerical  treat¬ 
ments  presented  and  will  be  discussed  for  the  granular  aggregate  case. 

Dusty  Gas 

For  a  dusty  gas  of  rigid  aluminum  particles  the  equilibrium  and  frozen 
(I)  sound  speeds  have  been  computed  (Ref  11).  The  frozen  (II)  sound  speed 
is  approximately  that  for  the  gas  alone.  The  numerical  results  are  shown 
In  Figure  3  for  densities  up  to  1 . 1  g/ cm^ .  The  frozen  (I)  and  equilibrium 
sound  speeds  were  close  enough  to  permit  plotting  only  one  curve  to  rep¬ 
resent  both.  The  frozen  (II)  sound  speed  is  about  330  m/s  by  definition. 
Higher  densities  were  not  considered  because  the  medium  at  1. 1  g/cm3 
is  already  beyond  the  dusty  gas  range,  as  shown  by  microphotographs  of 
the  powder.  The  density  of  the  medium  in  the  experiments  was  above 
1.3  g/cm^  showing  that  the  dusty  gas  approximation  is  of  limited  value. 

It  does  help  to  show  that  as  particles  are  added  to  a  gas  the  mass  loading 
causes  the  sound  speed  to  rapidly  decline,  flattening  out  at  a  low  value. 
The  sound  speed  only  increases  when  particle  contacts  sustain  load  and 
thus  increase  the  elasticity  of  the  medium. 

Porous  Elastic  Solid 


The  analog  to  the  frozen  (II)  sound  speed  in  the  dusty  gas  is  the  speed 
in  the  material  of  which  the  particles  are  composed.  Just  as  in  the  dusty 
gas,  there  is  some  small  quantity  of  energy  that  finds  its  way  through  a 
pure  gas  path,  so  for  the  porous  elastic  solid,  there  will  be  some  energy 
transmitted  through  a  purely  solid  path.  In  the  latter,  the  frequency  will 
have  to  be  high  enough  so  that  there  is  no  coupling  to  the  interstitial  gas. 
The  frozen  (H)  sound  velocity  for  an  aluminum  powder  is  therefore  about 
5180  m/s. 


The  equations  for  frozen  (I)  and  equilibrium  sound  speed  are  (Ref  12) 

,,  (2/^*  ^  ^ 

fo' 


■) ' (¥-) 


(7) 


(8) 


For  convenience,  it  has  been  assumed  that  the  interstitial  gas  behaves  like 
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a  perfect  gas  and  the  elastic  matrix  behaves  like  a  Hookian  material.  A 
single  prime  represents  the  solid,  a  double  prime  the  air.  Y  is  the  mass 
fraction  of  the  phase.  ^  is  the  initial  porosity  defined  as  the  fraction  of 
the  total  volume  occupied  by  the  gas.  at  ^  ^  specific  heat 

per  unit  volume  (at  constant  volume).  and ^  are  Lame's  constants  and 
^  is  the  thermoelastic  constant.  These  can  (Refs  13,  14,  15)  be  related 
to  common  mechanical  and  thermal  parameters  such  as  bulk  modulus 
and  coefficient  of  thermal  expansion.  Mumerical  values  are  available 
(Ref  16). 

Equation  7  states  that  the  square  of  the  frozen  sound  speed  in  the  mix¬ 
ture  is  exactly  equal  to  the  mass  fraction  weighted  mean  of  the  square  of 
the  isentropic  sound  speeds  in  the  elastic  matrix  and  the  interstitial  gas. 
For  the  mixture  Y'  =  0. 9995974  and  Y"  =  0.  0004026  and  tabulated  values 
for  the  sound  speed  in  aluminum  and  air  are  5182  m/s  and  331  m/s.  It 
follows  that  the  frozen  sound  speed  is  about  that  of  the  solid  of  which  the 
particles  are  composed.  Thus  both  frozen  sound  speeds  are  about  5180 
m/s  for  an  aluminum  powder.  For  a  mixture  of  aluminum  and  potassium 
perchlorate  there  may  be  a  narrow  range  of  velocities  determined  by  the 
properties  of  the  two  materials  and  the  statistical  distribution  of  paths 
for  the  sound  energy. 

Equation  8  was  used  to  calculate  the  equilibrium  sound  speed  for 
densities  above  1.1  g/cm^.  The  results  (shown  in  Figure  3)  indicate 
that  for  densities  above  1. 1  the  sound  speed  is  approximately  constant  at 
6400  m/s.  The  numerical  value  obtained  is  associated  with  the  data  used. 
The  point  to  note  is  that  even  for  a  density  as  low  as  1. 1,  a  very  high 
equilibrium  sound  velocity  is  calculated  as  soon  as  one  uses  the  porous 
elastic  medium  approach. 

Loosely  Packed  Granular  Aggregate 

Since  the  dusty  gas  and  porous  elastic  solid  approaches  provide  only 
extreme  limits  to  the  true  nature  of  the  medium,  and  give  sound  speeds 
so  far  apart,  and  bracketing  the  observed  values,  it  becomes  necessary 
to  try  and  make  a  better  calculation. 

Let  (T^  be  the  force  per  unit  area  applied  to  a  piston  of  area  Ap  com¬ 
pressing  unit  mass  of  the  mixture  a  distance  dx  (dV  =  Apdx) ,  then 
(T  *  r'-t-<r"  where  <r'  <r"  are  the  forces  per  unit  area  of  mixture  sup¬ 
ported  by  the  granules  and  gas  respectively.  The  equilibrium  soimd 
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where  Sq  is  the  total  entropy  at  T  =  Tq,  V  =  Vq  specific  volume  for  the 
mixture. 


If  Vq,  Vq'  are  the  initial  volumes  per  unit  mass  of  the  mixture  occupied 
by  the  mixture  and  granules  respectively,  at  temperature  Tq,  then 


where  F  is  a  function  to  be  determined.  Write  Fq  for  F  with  T  =  Tq. 

Now  Fq  has  been  measured  (Ref  1)  by  compressing  aluminum  powder  so 
slowly  that  the  interstitial  gas  remained  in  pressure  equilibrium  with  the 
atmosphere  and  an  isothermal  condition  was  achieved  at  the  temperature, 
Tq,  of  the  compression  device,  for  the  mixture.  The  results  are  plotted 
in  Reference  1  and  given  here  in  Table  1  together  with  least  square  fits 
for  1.4. 


The  role  of  temperature  in  equation  lO  is  assumed  to  be  primarily  that 
of  changing  the  size  of  the  granules  and  so  increasing  the  stress  with 
temperature.  This  assumption  is  reasonable  for  densities  above  1.4  g/cm3. 
An  increment  dT  will  increase  the  stress  by  dif' the  amount  required  to 
recompress  the  mixture  to  the  volume  that  existed  prior  to  raising  the 
temperature. 


dr' 


<i(T-T.)  =  - 


(11) 


Where  o.'dLT"  ^'wlth  a'  the  coefficient  of  thermal  e;q>ansion  of  aluminum 
and  zj  z'  the  specific  volumes  of  aluminum  at  To .  T  respectively. 

In  deriving  the  above  dCf'^O  and  dL^"  Ca  O  have  been  used  as 
characterizing  a  granule  thermal  expansion  with  restoration  to  same  mix¬ 
ture  volume  by  increase  in  applied  load.  Using  a  constant  coefficient  of 
thermal  expansion  a'  ,  and  integrating  from  Tq  to  T,  one  finds 

Vi  /  (12) 
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(13) 


For  (T  in  equation  9  we  can  use  the  perfect  gas  law 

where  V  is  the  volume  occupied  by  the  granules. 

It  can  be  shown  (Ref  17)  that  the  entropy  of  the  mixture  is 


where 


U(T) 


If  in  equation  15  we  use  U  -  cj'  T  and  define 

s  Y*c^  Y*'  C^"  so  that  U  ~  and  substitute  in 


and  substitute  in 


equation  14,  we  find  the  condition  for  S  =  Sq  is 

T  -  ,  o 

X>  ^ 

Equation  16  is  used  in  equations  12  and  l3.  These  are  then  differentiated 
and  used  in  equation  9  to  find  the  equilibrium  sound  speed.  However, 
equations  13  and  16both  contain  V.  It  is  assumed  that  V  fSg  V^'  where 
Vq'  is  the  initial  volume  occupied  by  the  granules.  The  justification  is 
that  the  principal  mechanism  in  the  deformation  range  considered  is 
rearrangement  of  granules  from  one  load-sustaining  structure  to-  another. 
The  density  range  for  this  to  be  valid  is  estimated  as  ^  p  f  /.  ? 

The  final  result  for  the  equilibrium  sound  speed  is  then 

rz'R 

'  li*  /  V.  LVi^) 

+  .^1  +  \i.  I'?  * ') 


The  results  of  calculation  with  this  formula  are  given  in  Table  2  and 
plotted  in  Figure  3.  As  the  density  is  lowered  from  1. 45  the  sound  speed 
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drops  rapidly  toward  the  dusty  gas  range.  This  drop  represents  a  de¬ 
crease  in  the  number  of  load-bearing  granule  contacts.  As  the  density 
increases  beyond  ^  =  1.6  the  sound  speed  rises  slowly  as  granular 
structure  rearrangements  occur.  One  expects  that  at  some  density  near 
that  of  aluminum  ( ^  =  2.7)  the  structure  will  become  locked  and  the 
sound  speed  will  simultaneously  rise  rapidly  to  that  of  the  porous  elastic 
solid  model. 

Discussion  of  Sound  Speed  Results 

In  the  introduction  to  this  Sound  Speed  Considerations  section,  it  was 
observed  that  to  determine  whether  the  observed  combustion  propagated 
Bubsonically  or  supersoncially  one  required  sound  speed  data.  Figure  3 
summarizes  the  results  of  the  calculations  and  also  shows  the  range  of 
speeds  observed  for  propagation  of  combustion  in  the  mixture.  The 
latter  are  obviously  supersonic,  leading  to  the  conclusion  that  the  observed 
propagation  of  combustion  is  not  a  simple  deflagration.  For  a  deflagra¬ 
tion  preceded  by  a  precompression  shock,  the  deflagration  would  have  to 
be  subsonic  with  respect  to  the  precompressed  medium.  Suppose  a 
density  change  from  l.  5  to  i.  8  across  the  shock,  then  by  conservation  of 
mass  the  compressed  mixture  directly  behind  the  shock  will  be  moving 

1  5 

in  the  same  direction  at  (l  -  7^)  times  the  shock  velocity.  For  a  shock 

at  900  m/s,  the  medium  will  be  traveling  at  750  m/s.  The  deflagration 
would  have  to  travel  at  900  m/s  relative  to  the  observer  to  maintain  a 
stable  profile.  However,  this  requires  the  deflagration  to  travel  at  750 
m/s  relative  to  the  compressed  medium.  From  the  sound  speed  results 
at  ^  =  1.8,  this  is  not  subsonic  as  required.  Therefore,  the  higher 
speed,  elongated,  stable  profile  can  not  be  a  deflagration  with  a  precom¬ 
pression  shock.  This  leads  to  the  conclusion  that  all  the  observed  pro¬ 
pagation  of  combustion  is  detonative,  which  is  the  same  as  the  conclusion 
reached  by  means  of  the  diffusivlty  considerations.  There  were  sufficient 
assumptions  (e.  g.  metal  dust  representing  the  mixture)  made  in  the  sound 
speed  calculations  to  raise  some  doubt  as  to  the  confidence  to  be  placed 
solely  in  the  absolute  values  of  the  sound  speed  results.  However,  they 
are  only  being  used  as  one  indicator  in  the  analysis  of  the  observed 
phenomena.  In  the  next  section,  still  another  ^proach  is  used  to 
examine  the  nature  of  the  observed  propagation  of  combustion. 
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ONE-PARAMETER  HYDRODYNAMIC  COMBUSTION  THEORY 


A  direct  approach  to  determining  whether  the  experimental  results 
represent  a  deflagration  is  to  calculate  the  speed  that  a  deflagpration  front 
would  have  in  this  particular  mixture.  However,  available  theories 
(Refs  29,  30)  must  assume  a  reaction  rate  law  to  provide  a  complete  set 
of  equations.  (For  a  detonation,  the  Chapman-Jouguet  stability  condition 
is  available.)  Introducing  this  assumption,  which  limits  the  validity  of 
the  results  (Ref  30),  is  not  necessary  for  the  purpose  at  hand.  If  the 
entire  deflagration  branch  of  the  Hugoniot  could  be  shown  to  lead  to  no 
propagation  velocity  comparable  with  that  observed  e:q>erimentally,  then 
it  would  have  been  proven  that  the  observed  combustion  was  not  a  deflag¬ 
ration.  The  one-parameter  hydrodynamic  combustion  theory  (OPHCT) 
which  is  developed  in  this  report  calculates  the  deflagration  branch  of  the 
Hugoniot  using  the  downstream  temperature  as  a  parameter.  Other 
features  of  the  Hugoniot  such  as  the  dead  zone  and  a  detonation  branch  of 
limited  validity  are  also  developed.  In  agreement  with  the  diffusivity  and 
sound  speed  results,  it  is  found  that  the  experimental  observations  can 
not  be  explained  as  a  deflagration.  In  the  remainder  of  this  section  the 
OPHCT  is  developed.  The  next  section  then  considers  the  detonation 
possibility. 

Basic  Assumptions  and  Equations 

We  shall  immediately  assume  the  following; 

a.  All  upstream  constituents  are  consumed  in  the  chemical  reaction; 
that  is,  combustion  is  complete. 

b.  No  external  heat  losses  or  external  work  need  be  considered; 
that  is,  the  problem  is  one  dimensional. 

c.  A  steady  state  exists  as  viewed  by  an  observer  traveling  with 
the  flame  front;  that  is,  there  is  a  constant  reaction  zone  profile  traveling 
at  constant  speed. 

d.  The  decomposition  and  interaction  of  K  Cl  can  be  neglected. 

e.  The  temperature  up  to  the  Chapman-Jouguet  plane  is  above 
2800**  K  so  that  it  can  be  assumed  that  liquid  A1  is  not  present. 

f.  The  species  chosen  as  products  are  the  significant  ones. 
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Al202«  A1  CI3,  and  Al203(g)  are  omitted. 

g.  An  appropriate  equation  of  state  is  that  for  a  perfect  gas  with 
the  volume  occupied  by  condensed  phase  AI2O3  subtracted  out.  The  equa¬ 
tion  of  state  of  condensed  AI2O3  is  ignored. 

Assumption  a  avoids  the  possibility  that  some  of  the  aluminum  is  not 
consumed,  which  would  introduce  a  second  parameter  and  require  con¬ 
sideration  of  the  thermal  energy  sink  action  (Ref  31)  of  the  unconsumed 
particles.  Since  full  consumption  of  the  aluminum  will  lead  to  the  highest 
propagation  rate,  assumption  a  was  justified. 

Assumptions  b  and  c  make  the  problem  an  algebraic  one  by  permitting 
use  of  the  discontinuous  solutions  (jump  conditions)  of  the  differential 
equations  representing  conservation  of  mass,  momentum,  and  energy. 
Appendix  A  derives  the  jump  conditions  in  a  manner  which  shows  that 
assumption  c  is  needed  for  a  finite  thickness  flame  zone. 

Assumptions  d,  e,  and  f  are  of  the  type  one  must  always  make  in 
combustion  calculations;  i.e.  one  must  choose  the  products.  These 
choices  were  guided  by  experience  in  the  aluminized  propellant  field 
since  the  primary  Interest  was  in  deflagrations  (low  pressures).  In  fact, 
these  choices  serve  to  invalidate  the  calculation  of  the  detonation  branch 
by  the  OPHCT  because,  as  will  be  evident  in  the  detonation  section  of  this 
report,  other  products  must  be  assumed  for  the  high  pressures  of 
detonation. 

Assumption  g  is  again  valid  for  the  low  pressures  of  deflagrations 
and  invalid  for  the  detonation  branch  where  one  must  use  an  equation  of 
state  suitable  for  high  pressures  (Ref  32) . 

Next  we  state  the  applicable  equations. 

Let  U  be  the  velocity,  positive  to  the  right  of  the  combustion  zone. 

Then  if  Uj^  and  U2  are  the  particle  velocities  upstream  and  downstream 
within  the  stable  reaction  zone  profile,  and  v^  and  V2  are  the  same  rela¬ 
tive  to  the  combustion  zone,  we  have  v^  =  ui  -  U  =  -u  since  u.  =  0  and 
V2  =  U2  -  U. 

Following  a  mass  element  in  the  flow,  the  jump  conditions  (Rankine- 
Hugoniot)  as  derived  in  Appendix  A  are,  for  the  conservation  of  mass 


for  conservation  of  momentum 


for  the  conservation  of  energy 


(19) 


(20) 


These  can  be  combined  into  the  Bigoniot  Adiabatic  Equation 

=  '2  )  (21) 

which  replaces  the  energy  equation. 

Here  h,  and  functions  of  both  the  temperature  and  composi¬ 

tion,  and  under  certain  assumptions  could  also  depend  on  the  pressure. 


A  fourth  equation  that  applies  is  the  equation  of  state,  which,  if  one 
considers  the  volume  occupied  by  the  condensed  AI2O3,  becomes  the 
following  modification  of  the  perfect  gas  law. 

JL  -  ^  ^  ^ 

h  (22) 


^  - 


The  indicates  a  summation  over  all  species,  whereas  Z*  is  used 
for  gaseous  species  only.  The  and  yr%.  are  the  mole  fraction  and 
molecular  weight  respectively.  Using  the  same  notation,  the  enthalpy  of 
the  product  mixture  used  in  equation  2l  can  be  written  as 

r-x.  ^  C,,. 

s'  yt,rr,^  (23) 

Here  Cfi^  and  ^ao^are  molar  specific  heat  and  enthalpy  at  Tj  =  298. 15*  K 
respectively  for  each  product  species. 

The  mixture  is  60/40,  by  weight,  of  potassium  perchlorate/aluminum. 
The  loading  density  of  the  powder  is  1.  5  grams/cm  (  p,  =1.5).  Thus, 
considering  a  volume  of  1  cm^,  the  combustion  involves  the  following 
moles  (  and  )  of  reactants: 


L 


O.t’^hS 


OJL 
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The  chemical  reaction  that  will  be  used  is 


)n  -t-  Y>  n  -f  w 

ae.  Q/.o  cLf^o  0 


w 


All  products  except  AI2O3  (condensed)  are  considered  to  be  in  the  gaseous 
state  in  the  range  of  product  temperatures  of  interest.  For  forming 
mole  fractions  from  the  it  is  found  convenient  to  use  f/  as  defined 
below,  omitting  the  condensed  species,  as  the  normalizing  factor  and 
make  adjustments  as  required  for  the  condensed  AI2O3. 


(24) 

(24a) 


The  combustion  provides  three  material  balance  equations  (e.g. 
aluminum  can  only  go  into  aluminum-containing  species) .  These  are 


"at  *  "oho  *  ^”at^o  =  o.  ozii  ps) 


=  "*«,v  = 


(27) 


where  the  right  sides  are  the  number  of  moles  per  cm^  provided  initially. 
These  can  be  rewritten  in  terms  of  mole  fractions. 


X  X 
OJ. 
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(28) 

(29) 

(30) 


Since  we  consider  that  equilibrium  exists  among  the  products,  for 
each  molecular  species  involving  A1  and  O2  we  can  write  an  equilibrium 
reaction  and  an  equation  for  the  equilibrium  constant,  which  is  a  tabulated 
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function  of  temperature.  For  example: 


oe  ■'•iq 


cao 


hho  ■  s  k 


CUtO 


(31) 


'/t/  •  '^2. 

Note  that  the  use  of  O2  as  the  source  of  oxygen  results  in  fractional 
powers.  However,  the  available  tables  give  the  equilibrium  constants 
for  O2  as  the  oxygen  source.  We  can  circumvent  this  difficulty  as 
follows:  Consider  the  reaction  written  as  CUt  "*■  O  ^  OJiO  for  which 


oex> 


'Ota 


since 


OiO 


mo 


Thus  we  can  use  reactions  based  on  O,  eliminating  fractional  powers, 
except  for  O ,  by  combining  tabulated  eqpiilibrium  constants  in  ^prop- 
riate  fashion.  Further,  we  can  use  the  fact  that  the  partial  pressures 
s  to  write  the  equilibria  equations  in  terms  of  the  mole 

fractions.  The  vapor  pressure  of  Al203(c)  is  considered  negligible,  and 
is  included  in  the  tabulated  equilibrium  constant  for  Al203(c).  The 
resulting  equations  are: 


K  . 

aio 


OtO 


S  A- 


(32) 


X 

*  f 

OM.  O 
I 

"at  o 


•'o’ 


(33) 


(34) 


(35) 


The  functions  a,  b,  c,  and  d  are  defined  in  the  above  equations. 

Equation  2l  combined  with  23,  equation  24  or  24a,  and  equations  22, 
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28,  29,  30,  32, 
unknowns 


33,  34,  and  35  constitute  ten  equations  in  the  eleven 

XX  ^  ^  XX 

a  OfO  ktt  Oij)  O 


These  equations  represent  simple  conservation  of  mass,  momentum, 
energy,  equation  of  state,  material  balances,  and  equilibrium  conditions. 
There  is  one  more  unknown  than  equations  and  as  stated  the  approach  is 
to  use  T2  as  a  parameter  and  seek  a  solution.  It  was  found  possible  to 
reduce  this  system  of  equations  to  a  single  polynomial  equation  in 

ffe.’b)'  Since  the  reduction  procedure  with  appropriate  modi- 
fications^s  applicable  to  other  problems, it  is  given  in  detail  in  Appendix 
B.  For  each  temperature,  the  positive  real  roots  of  the  polynomial  which 
lead  to  positive  pressures  were  the  only  physically  significant  ones.  The 
problem  was  programmed  for  the  IBM  709  computer  and  after  these  roots 
were  found,  the  downstream  values  of  the  unknowns  were  calculated.  The 
values  obtained  were  unique  and  are  listed  in  Table  3  as  a  function  of 
temperature. 


Numerical  Results  and  Discussion 


The  values  in  Table  3  were  selected  from  the  full  computation,  which 
was  done  for  temperatures  100*K  apart.  As  expected,  one  finds  a  de¬ 
flagration  branch,  a  physically  impossible  dead  zone  (contradicts  con¬ 
servation  of  mass  and  momentum),  and  a  detonation  branch.  The  deton¬ 
ation  results  must  be  regarded  as  of  limited  validity  because  of  the 
assumptions  d,  e,  f,  and  g  made  at  the  beginning  of  this  section.  Three 
conclusions  follow  immediately. 

a.  Deflagration  can  occur  only,  for  temperatures  at  the  rear  of 
the  stable  reaction  zone,  below  4000* K  and  only  with  flame  speeds  of 
the  order  of  cm/sec. 

b.  The  dead  zone  extends  from  4000*  K  upwards  to  a  temperature 
approaching  7000*  K,  the  exact  value  depending  on  the  limiting  assump¬ 
tions  d,  e,  f,  and  g. 

c.  A  Chapman-Jouguet  detonation  can  be  e:q>ected  to  have  a  tem¬ 
perature  in  the  neighborhood  of  7000*  K. 

The  OPHCT  has  served  to  confirm  once  again  that  the  e;q;)erimentally 
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observed  phenomena  are  not  those  of  a  deflagration.  There  are  other 
extensions  of  OPHCT  (covolumes,  other  product  species,  incomplete 
combustion  of  Al)  which  could  be  carried  out  but  they  would  not  change 
this  conclusion.  For  this  reason,  we  turn  now  to  calculating  the  Chap- 
man-Jouguet  detonation  for  the  mixture, 

CHAPMAN-JOUGUET  DETONATION  CONSIDERATIONS 

The  calculation  of  the  speed  of  an  ideal  detonation  differs  from  that 
of  a  deflagration  in  that  the  system  of  equations  is  complete.  This  is 
achieved  by  the  addition  of  the  Chapman-Jouguet  stability  condition.  The 
problems  with  respect  to  the  detonation  calculation  stem  from  the  high 
pressures  (e.g.  100,000  atmospheres)  in  the  detonation  zone.  The  equa¬ 
tion  of  state  of  the  gaseous  products,  that  for  solid  products,  and  the 
choice  of  products  and  of  thermodynamic  data  must  all  be  appropriate  to 
the  anticipated  higher  pressures.  The  Chapman-Jouguet  values  are  cal¬ 
culated  by  computer.  At  Los  Alamos  Scientific  Laboratories  a  Brinkley- 
Kistiakowsky-Wilson  (BKW)  technique  is  used.  At  Lawrence  Radiation 
Laboratory  a  similar  program  which  can  handle  two  condensed  product 
species  is  used.  The  latter,  called  RUBY  (Ref  33)  was  made  available 
to  Picatinny  Arsenal  where  it  was  used  for  the  calculations  described  in 
this  report.  The  RUBY  program  uses  a  Kistiakowsky- Wilson  equation  of 
state  for  the  product  gases  (Refs  34,  47).  It  requires  one  to  choose  the 
products,  provide  thermodynamic  data  for  them  and  equations  of  state 
for  the  condensed  products.  The  choices  made  and  parametric  studies  of 
the  significance  of  these  decisions  are  given  below.  This  is  followed  by 
a  presentation  and  discussion  of  the  results. 

Inputs  and  Parametric  Studies 

The  products  chosen  differ  from  those  of  OPHCT  in  two  important 
regards.  The  gaseous  product  AI2O3  is  considered  to  exist  at  the  high 
pressures.  It  is  recognized  that  the  literature  (Refs  35-39)  indicates 
that  AI2O3  gas  decomposes,  producing  suboxides.  However,  although 
this  is  true  at  very  low  pressures  compared  to  those  of  a  detonation, 
experience  with  calculations  on  explosives  containing  aluminum  indicates 
a  need  for  AI2O3  gas  to  obtain  correct  results.  The  other  difference 
from  OPHCT  is  that  the  KCl  is  permitted  to  decompose  completely  and 
react,  thus  allowing  for  the  additional  species  K,  K2,  and  AICI3,  The 
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suboxides  of  aluminum  have  again  been  included.  The  condensed  species 
expected  are  A1  (s)  and  AI2O3  (s).  The  species  used  are  listed  in  Table 

5. 


Thermodynamic  tables  for  the  products  are  available  from  many 
sources  (Refs  40-46).  Here,  too,  one  must  be  careful  to  select  data 
appropriate  for  high  pressures.  The  RUBY  program  minimizes  the  free 
energy  to  find  the  chemical  equilibrium  among  the  products  for  the  T  and 
p  under  consideration.  The  free  energy  is  computed  from  specific  heat 
vs  tempeiature  tables.  If  such  tables  have  simply  assumed  a  constant 
specific  heat  for  high  temperatures,  this  will  affect  the  computation. 
Unfortunately,  this  occurs  occasionally  in  the  JANAF  tables  (Ref  43)  and 
for  such  species  the  JANAF  tables  were  not  used.  The  JANAF  data  were 
used  in  the  form  of  the  generating  pol3momials  (Ref  44) .  However,  most 
of  the  thermodynamic  data  were  taken  from  the  LASL  compilation  (Ref  45), 
which  has  a  history  of  successful  use  in  detonation  calculations.  The 
thermodynamic  data  used  in  the  calculation  are  listed  in  Table  5. 

The  covolumes  for  the  gaseous  species  are  (Refs  34,  47)  equal  to 
10,46  Vj,  where  Vj  is  the  molecular  volume  in  A^.  For  KCl  the  molecular 
volume  was  calculated,  assuming  full  rotatibility,  using  (Ref  48)  an 
equilibrium  internuclear  distance  of  2. 671  A,  and  Van  derWaal  radii  for  Cl 
and  K  of  1, 80  A  and  2, 825  A,  respectively.  The  last  value  was  calculated 
by  assuming  that  the  Van  derWaal  radii  can  be  taken  as  equal  to  covalent 
radius  plus  0, 80  A.  The  covolumes  used  are  shown  in  Table  5. 

The  equations  of  state  at  high  pressures  for  A1  (s)  and  AI2O3  (s)  were 
of  the  form  used  by  Cowan  and  Fickett  (Ref  34) . 

I>  *  a(V,)  r  »  k(Vs)r^ 

Here  Y  =  initial  specific  volume/  specific  volume  at  T  and  p.  The  Aj 
coefficients  used  were  those  of  the  shock  Hugoniots  for  A1  (s)  (Ref  49) 
and  AI2O3  (s)  (Ref  50)  used  as  an  adequate  approximation  to  this  part  of 
the  equation  of  state.  The  coefficients  Bq  and  B.  were  determined  using 
thermal  expansion  (Ref  51)  and  compressibility  (Uta  (Ref  52).  (The 
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former  data  were  ‘also  used  for  generating  a  thermal  expansion  ec[uation 
which  is  also  required.)  The  coefficients  Cq,  Cj^,  and  C2  are  associated 
with  electronic  vibrations  and  (for  the  pressures  obtained  with  this  mix¬ 
ture)  can  be  set  equal  to  zero.  The  coefficients  for  the  equations  of  state 
and  thermal  expansion  equation  are  listed  in  Table  5, 

The  principal  calculation  made  use  of  the  best  data  available.  How¬ 
ever,  other  runs  were  made  to  determine  how  the  results  depended  on 
the  data  used.  The  equations  of  state  were  deleted  and  the  solids  treated 
as  incompressible.  In  a  second  variation,  some  of  the  thermod3niamic 
data  from  LASL  were  replaced  by  JANAF  data.  The  effect  of  loading 
density  of  the  mixture  was  evaluated  by  changing  the  value  from  1.  5  to 
0.  5  g/cm3.  Dr.  C.  Mader  of  LASL  was  kind  enough  to  calculate  the  C-J 
values  with  his  BKW  code.  Since  several  less  important  species  were 
omitted  by  him,  this  provided  both  a  measure  of  their  effect  and  by  virtue 
of  the  agreement  with  the  principal  RUBY  calculation  a  corroboration 
thereof. 

Numerical  Results 

The  results  of  the  RUBY  calculations  of  the  Chapman-Jouguet  detona¬ 
tion  for  the  mixture  are  shown  in  Table  6.  Although  a  variety  of  input 
conditions  are  represented,  the  results  all  show  that  the  calculated 
detonation  velocity  is  several  times  the  highest  experimentally  observed 
value. 

From  the  results  one  can  also  conclude  the  following;  The  equations 
of  state  should  be  included  in  any  calculation.  Both  A1  (s)  and  Al'203  (s) 
are  required.  Although  some  gaseous  species  (e.g. ,  AI2O2)  could  be 
omitted  for  this  mixture  and  others  such  as  AI2O3  (g),  Alcfs  (g)  should 
not  be  omitted,  it  would  be  best  for  a  similar  mixture  to  include  all  those 
listed.  The  thermodynamic  data,  as  eiqpected,  affect  the  results.  A 
reduction  in  density  reduces  the  detonation  velocity  and  pressure  and 
AI2O3  (s)  is  no  longer  present. 
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OVERALL  DISCUSSION 


The  sections  on  diffusivity  and  sound  speed  considerations  and  the 
OPHCT  theory  established  that  the  propagation  of  combustion  was  not 
deflagrative.  The  calculated  Chapman-Jouguet  detonation  velocity  was 
about  3  to  4  times  the  highest  experimentally  observed  speed  of  the  flame 
front.  However,  this  last  result  only  indicates  that  we  are  not  dealing 
with  an  ideal  detonation. 

In  the  ideal  detonation,  the  problem  is  considered  as  strictly  one 
dimensional  with  a  reaction  zone  of  negligible  thickness.  The  experi¬ 
mental  results  showed  a  reaction  zone  several  centimeters  long  and  con¬ 
siderable  lateral  rarefaction  starting  within  one  centimeter  of  the  front. 
Hence  one  possible  explanation  of  the  discrepancy  between  the  RUBY  and 
eiqierlmental  results  is  the  lateral  rarefaction. 

The  RUBY  calculation  is  not  designed  for  calculations  that  include 
relaxation  phenomena  associated  with  physical  kinetics  of  the  mixture. 

In  the  mixture,  there  are  aluminum  particles  that  require  times  to 
vaporize  that  are  significant  fractions  of  the  time  required  by  the  reaction 
zone  to  pass  a  fixed  point.  In  the  experimental  results,  this  slower 
release  of  the  fuel  would  affect  the  flame  speed.  This  is  in  accord  with 
the  previously  reported  (Ref  4)  dependence  of  combustion  on  particle  size, 
and  the  observed  structure  of  the  reaction  zone.  The  presence  of  two 
reaction  zone  profiles  with  different  speed  ranges  also  Indicates  the 
presence  of  some  rate  phenomena.  Similar  physical  effects  have  been 
reported  for  nitrocellulose  or  nitroglycerin  powder  (Ref  53) ,  liquid 
nitroglycerin  or  methyl  nitrate  (Refs  54,  55),  flaked  TNT  and  composi¬ 
tion  B  (Refs  56,  57  ,  58),  and  others  (Ref  59).  Hence  another  possible 
explanation  of  the  divergence  of  the  RUBY  results  from  the  observations 
is  that  one  must  include  physical  kinetics  in  the  calculation. 

Further  Experimental  Work 

The  effect  of  lateral  rarefaction  could  be  evaluated  by  altering  the 
confinement  of  the  column  of  mixture  and  using  larger  diameter  columns. 
Other  worthwhile  modifications  include  using  longer  propagation  sticks, 
var3ring  the  mode  of  Initiation,  and  examining  the  reaction  zone  by  time- 
resolved  spectroscopy.  In  addition,  it  would  be  worthwhile  reviewing 
the  experimental  woiic  on  low  spaed  detonations  to  observe  the  techniques 
used  by  others. 
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Further  Theoretical  Work 


Using  the  experimental  data  to  give  a  measure  of  the  lateral  rare¬ 
faction  rate  and  ignoring  the  relaxation  effects  of  the  medium,  one  can 
calculate  estimates  of  the  lateral  effect  on  the  detonation  velocity.  This 
could  then  be  compared  with  the  experimental  results. 

The  area  in  which  it  is  believed  there  is  a  greater  contribution  to  be 
made  is  that  of  the  relaxation  effects.  Prior  to  working  on  this  problem, 
it  would  be  necessary  to  review  in  detail  studies  of  reaction  zone  structure 
(Refs  60-64)  with  a  view  to  substituting  or  adding  physical  kinetic  consid¬ 
erations  (vaporization,  diffusion)  to  the  chemical  reaction  rates  considered. 
Another  approach  would  be  based  on  a  two-phase  flow  theory  (Ref  65) 
where  the  solution  would  depend  on  the  assumptions  made  on  mass  and 
energy  interchange  and  forces  acting  between  the  two  phases.  Finally, 
one  can  explore  the  concepts  used  for  detonations  in  dilute  sprays  (Ref 
66).  An  explanation  of  low  speed  detonations  and  the  experimental  obser¬ 
vations  on  this  mixture  would  be  the  objectives  of  the  theoretical  work. 
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Table  1 


Data  on  the  Compression  of  Aluminum  Powder*** 


V  -  V** 

Fo* 

V  -  V** 

Fo* 

— V — 

''o 

0. 

0. 

0.1178 

75.  01 

0.0036 

0.01379 

0. 1206 

78.60 

0.0078 

0.09652 

0.1228 

82. 18 

0.0120 

0. 2620 

0. 1261 

85.  21 

0.0162 

0. 5515 

0.1283 

89. 19 

0. 0196 

0.9928 

0. 1293 

92.93 

0.0230 

1.599 

0.1331 

96.80 

0.0271 

2.303 

0.1363 

100.10 

0. 0298 

3.185 

0.1395 

104.78 

0.0332 

4. 247 

0.1422 

107.83 

0.  0398 

6.  673 

0. 1469 

115.  83 

0.0424 

8.053 

0. 1495 

121.07 

0. 0463 

9.652 

0.1520 

124.1 

0.0528 

12.961 

0. 1556 

129.06 

0. 0541 

14.89 

0.1581 

132.37 

0. 0579 

16.82 

0.1602 

137.34 

0.0610 

18.75 

0.1632 

140.65 

0.0635 

20.96 

0.1657 

144. 51 

0, 0673 

23.33 

0.1682 

148. 92 

0.0691 

24.93 

0.1706 

153.88 

0.0759 

30.23 

0.1736 

157.19 

0.0789 

32.96 

0.1755 

162.43 

0.0820 

35. 58 

0.1789 

165.47 

0.  0850 

37.89 

0.1808 

170.98 

0.0879 

40.54 

0.1832 

175.12 

0.0903 

43.30 

0.1859 

179.25 

0.0933 

46. 05 

0.1880 

184.77 

0.0962 

49.09 

0.1905 

189.18 

0.0997 

52.12 

0.1929 

193.04 

0.1049 

58.18 

0.1950 

198.56 

0.1077 

61.50 

0.1975 

203.52 

0.1100 

65. 08 

0.1999 

207.38 

0.1128 

68.11 

0. 2024 

212.35 

0.1150 

71.49 

0.2045 

217.86 
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Table  1  (Continued) 


Data  on  the  Compression  of  Aluminum  Powder’*""* 

Vn  -  V** 

Fo* 

Vo  -  V"-" 

Fo* 

Vo 

Vo 

222. 00 

0.2074 

0.2205 

252.34 

226.14 

0. 2094 

0. 2234 

256. 47 

231.65 

0.2119 

0.  2251 

261. 99 

237.17 

0.2135 

0.  2271 

267. 50 

242.  68 

0.2156 

0. 2300 

271.64 

247.  65 

0.2181 

*  Values  of  Fg  (pressure)  shoum,  after  multiplication  by  10^,  are  in 
c.g.s.  units. 

**  A  reference  of  Vq  =  l/^  =  been  used  as  the  starting  point  for 

significant  pressures.  The  actual  pressure  was  0. 0120  dynes/cm^. 

***  Least  square  fits  using  polynomials  up  to  and  including  the  7th  degree 
were  tried  for  the  entire  data  and  with  the  data  divided  in  two  parts.  The 
results  chosen  as  most  suitable  for  sound  speed  calculations  were 

for  X  r  ^  0.05786  Pq  =  Aj,+Aj^X+A2x2+A3X3+A4X^+A5X® 

with  Aq  =  -0.  871223  x  10“2 

A^  =  -0.311934  X  10^ 

A2  =  0. 194468  X  102 

Ag  =  0.330135  x  10® 

A4=  -0.524995  X  10“^ 

A5  =  0.325772  x  10® 


for  Xa.^°  ~  ^  >  0.05786 
Vo 

with  Aq  =  0. 234946  x  102 
Ai  =  *0.  942820  x  10® 
A2  =  0. 178666  x  10® 


Fq  =  Ao+AjX+A2x2+A3X®+A4X^ 

A„  =  -0.  658572  x lO® 

A4  =  0. 115540  X  10® 
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Table  2 

Sound  Speed  Restdts  for  the  Granular  Aggregate  Model 


Density  (g/cm®) 


l^ed  (m/sec) 


1.40 

37. 

1.45 

151 

1.50 

206 

1.55 

250 

1.60 

269 

1.65 

280 

1.70 

289 

1.75 

299 

1.80 

312 

1.85 

328 

1.90 

346 

1.95 

365 

2.00 

386 

2.05 

407 

2.10 

427 

2.15 

447 

2.20 

466 

2.25 

484 

2.30 

499 

2.35 

513 
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Table  3 


One-Parameter  Hydrodynamic  Combustion  Theory  Results 


T*K 

U  (cm/s) 

V2  (cm/s)  ^2  (gm/cm2) 

Deflagration  Branch 

/^2  (bars)* 

3000 

2.26 

2.92  X  10® 

1.16  X 10-5 

0.02 

3500 

3.70 

1.  57  X  105 

3.  53  X  10-5 

0.14 

3900 

0.61 

4.  94  X  103 

1.84  X 10-4 

1.00 

Dead  Zone  (i>9  ^  but  ^  9  <  (1*5) 

4000 

8.48 

4.  57  X  104 

2.  79  X  10-4 

1.57 

7000 

2.45  X  10® 

2.  68  X  10®  1.  38 

Detonation  Branch 

8.  09  X  10® 

7100 

3.8  X  10® 

3.63  X  105 

1.57 

9.  52  X  10® 

7600 

1.76  X  10® 

1.05  X  10® 

2.51 

1.86x104 

TK 

^A1 

^AlO  ^Al20 

^2 

3000 

2.0  X  10-2 

1.8  X  10-1 

3.2x10-2  1.9x10-1 

6.8  X  10-4 

3500 

2.2  X  10-1 

1.  6  X  10‘1 

1.  lx  10-1  1.6x10-1 

2.9  X  10"® 

3900 

2.5  X  10"! 

1.5  X  10-1 

1.3xl0"l  1.7x10-1 

4. 4  X  10-2 

7100 

5.8  X  10-2 

7. 1  X  10-2 

1.2x10-1  2.3x10-1 

8.4x10-4 

7600 

3.7x10-2 

1.0  X  10*1 

1.5xl0"l  1.9xl0"l 

1.7  X 10-® 

^KCl  ^Al203(c) 


5. 

GO 

X 

10-1 

6.9 

X 

O 

1 

•-* 

CO 

1 

X 

10-1 

2.3 

X 

1 

o 

2. 

5 

X 

10"^ 

1.3 

X 

o 

1 

5. 

1 

X 

10"^ 

5.5 

X 

1 

o 

5. 

2 

X 

1 

O 

5.6 

X 

1 

o 

*  1  atmosphere  =  1. 01325  bars  »  1. 01325  x  10^  dynes/ cm^ 
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Table  4 


Tentative  Hugoniot  Data  for  AI2O3  (s)  * 


P  (meKabars) 

Y** 

f 

0 

1 

0.02655 

1.0141 

0.05528 

1.0275 

0.08619 

1.0401 

0.11927 

1.0522 

0. 15453 

1.  0637 

0.19193 

1.0746 

0.23157 

1.0850 

0.27339 

1.0948 

0.31761 

1.1042 

800285E-5  A^  =  -0.428736E2 

Ao  =  0. 131033E3 

Ag  = -0. 13509E3 

A^  =  0. 469329E2 

in  f5  =  Aq  +  A^Y  +  A2Y^  +  AgY^  + 
by  least  squares  fit 


*  These  preliminary  data  were  supplied  in  different  form  as  a  courtesy 
by  N.  Coburn,  Naval  Ordnance  Laboratory,  Silver  l^rings,  Md. 

Y  =  0.27027/  specific  volume  at  |>  of  shock. 
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Table  5 


RUBY  faput  Data 

Thermodynamic  Data 


Species 

Specific 

Heat 

Source 

Covolume* 

Heat  of 
Formation** 

Free 

Energy  of 
Formation** 

A1 

LASL 

350 

3.2659 

2. 8597 

AlO 

DOW 

1160 

0.8959 

0.  6350 

AlgO 

DOW 

1300 

-1.3164 

-1.6121 

AI2O2 

DOW 

1800 

-3.9943 

-4. 0258 

AI2O3 

LASL 

1350 

-10. 0769 

-9.8861 

AICI3 

LASL 

2600 

-5.8795 

-5.7366 

K 

LASL 

920 

0.8998 

0.6123 

K2 

LASL 

4600 

1.2869 

0. 9222 

K  Cl 

LASL 

2800 

-2.1605 

-2.3530 

O2 

DOW 

350 

0.0 

0.0 

0 

DOW 

120 

2.4937 

2.3194 

Al  203(8) 

LASL 

24.485 

-16. 7647 

-15. 8301 

Al(s) 

LASL 

9.9889 

0.0 

0.0 

*  In  cm^/mole.  For  solids,  value  is  volume  in  cm®  occupied  by  one 
mole  at  initial  T  and  jy  . 

**  In  megabar  cm®  psr  mole. 
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Table  5  (Continued 
Thermodynamic  Data  (Continued) 


Species 

GAM  0* 

GAM  1 

GAM  2 

GAM  3  GAM  5 

A1 

2.4135997 

-8.7089276E- 

-3  7.1777077E-3 

0.  0. 

AlO 

4.2483855 

1.1573301 

-1. 8338290 

1.5959543  -7.  0816294E-4 

AlgO 

6. 5471740 

1.5580869 

-4. 1013633 

3.  5731770  -1.  585961 6E-3 

AI2O2 

9.  2748740 

2.8285952 

-7. 4501969 

6.4935588  -3.  2106461E-3 

AI2O3 

12.324936 

0.71581680 

-0.6036318 

0.  0. 

AICI3 

9.  5651999 

0.2642911 

-0. 2215514 

0.  0. 

K 

2.  072784 

2.0914092 

-1.4225951 

0.  0. 

^2 

3.4132741 

8.4193720 

-3.942406 

0.  0. 

KCl 

4.3330734 

0.6990020 

-1.9300586E-3 

0.  0. 

O2 

3.8926835 

3.0720323 

0. 50698165 

-4.  6228122  -7. 805438E-4 

0 

2.4786804 

-0.91630437 

3.4551387 

-2.485581  3.7995423E-5 

Al203(s) 

13.77977 

3.965481 

-7.377259 

4. 357859  0. 

Al(s) 

2. 8406157 

0.32750082 

-0. 4891151 

0.  0. 

Q 

*  The  specific  heat  is  in  megabar  cm  /e.v.  mole  and  is  given  by  the 
following  expansion  with  T  in  electron  volts.  ^ 

=  GAM  0  +  GAM  1  T+GAM  2  T*+  GAM  3  T  +  GAM  4  T*  +  GAM  5  T  +  GAM  6 

The  GAMOS  coefficients  have  corresponding  units.  All  GAM  4  and  GAM  6 
were  zero  for  the  data  used.  (E2  means  x  10^,  etc.) 
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Table  5  (Continued) 


Equation  of  State  Data* 

^0 

^3 

^4 

Al  203(8) 

-0.  80028 5E -5 

-0. 428736E2 

0. 131033E3 

-0. 13509E3 

0. 469329E2 

Al(s) 

0.466 

■0.1269E1 

0.375 

0.428 

Bo 

®1 

ALPS 

BETS 

GAMS 

Al  203(0) 

0. 1405 

0.22687 

1. 0092533 

-0.36321 

0. 12721 

Al(s) 

-0. 16284 

0.39467 

1.0261949 

-1.  04433 

0.98718 

*  The  pressure  is  in  megabars  and  is  given  as  a  function  of  Y  =  initial 
volume/volume  at  T  and  p  . 

j>  =  Aq+  AjY  +  +  AjY^  +  A4Y^  +  Bq  +  B^Y 

The  thermal  equation  of  state  relates  Y  to  T  at  atmosphere  pressure  by 
Y  =  ALPS  +  BETS  •  T  +  GAMS  •  T* 
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Table  6 

Results  of  RUBY  Rons** 


Run  No.  55  56  47  48  LASL  45 

Eq.St.Al203(s)  *  INCOMP  INCOMP  INCOMP  INCOMP  INCOMP 
Eq.St.Al(B)  *  *  *  ^  __  inCOMP 


Thermo  data 


Species 

Density 

* 

* 

* 

* 

T*K 

6267 

8731 

^  megabars 

0.1070 

0.0944 

0. 

U  m/s 

4966 

4977 

Al 

3.9E-5 

1.2E-4 

7. 

AlO 

2.1E-7 

2.8E-6 

1. 

AI2O 

7.4E-5 

5.7E-5 

9. 

AI2O2 

l.lE-7 

2.1E-7 

2. 

AI2O3 

4.0E-3 

4. 4E-4 

5. 

AICI3 

1.3E-3 

8.3E-4 

4. 

K 

3.8E-3 

2. 4E<3 

1. 

Ko 

9.9E-7 

6.1E-5 

2. 

KCl 

5. 6E-4 

1.  9E-3 

3. 

O2 

2.1E-9 

7.1E-8 

2. 

0 

6.4E-8 

1.8E-6 

1. 

Al203(8) 

1.7E-3 

5.3E-3 

0. 

Al(s) 

1.9E-3 

2.3E-3 

8. 

*  DOW  *  DOW 


* 

0.5 

* 

* 

REDUCED  REDUCED 

*  * 

5775 

7372 

7290 

7372 

0167 

0.0761 

0. 1026 

0.  0407 

3115 

4545 

4845 

3286 

OE-4 

3.5E-7 

1.9E-3 

3.2E-6 

7E-5 

1.2E-7 

— 

3.2E-6 

7E-4 

7.5E-8 

~ 

7.4E-7 

lE-6 

3.6E-9 

— 

9.1E-8 

4E-3 

1.6E-4 

1.7E-3 

— 

lE-4 

5.2E-5 

1.3E-3 

~ 

2E-3 

1.5E>4 

4.2E-3 

lE-5 

2.1E-6 

1.8E-6 

-- 

lE-3 

4.  2E-3 

1.7E-4 

4.3E-3 

3E-8 

l.OE-6 

— 

1.2E-4 

3E-6 

2.  OE-6 

— 

9.8E-5 

5.  6E-3 

4.  lE-3 

5. 7E-3 

8E-4 

3.2E-3 

— 

3.5E-3 

Indicates  that  standard  conditions  as  described  in  report  and  Table  5 
were  used.  INCOMP  indicates  that  equation  of  state  was  deleted  and  the 
solid  considered  incompressible.  A  run  with  density  0. 5  g/cm^  instead 
of  the  usual  1. 5  was  made.  DOW  indicates  that  the  thermodynamic  data 
from  LASL  for  Al,  AICI3,  Al203(s),  Al(s)  was  replaced  by  that  of  DOW. 
REDUCED  indicates  that  the  species  indicated  by  —  in  the  data  colunrm 
were  deleted. 

is  the  detonation  velocity,  are  the  Chapman-Jouguet  pressure 
and  temperature,  ^cie  concentrations  are  per  gram  of  initial  mixture. 
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Table  6  (Continued) 

The  run  numbers  are  local  designations  and  are  convenient  for  reference 
and  making  comparisons.  LASL  is  the  run  made  at  Los  Alamos  Scientific 
Laboratory  with  the  BKW  code  by  Dr.  C.  Mader.  There  are  a  few  other 
differences  between  runs  which  are  not  indicated  in  the  tabulation  because 
the  effect  on  the  results  was  shown  separately  to  be  of  no  significance. 
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Fig  1  Propagation  of  Combustion  in  the  Granular  Mixture 


Fig  2  Two  Reaction  Zone  Profiles 
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(4/W>  p“"*$ 


Sound  Speed  Considerations 


Appendix  A 


DERIVATION  OF  JUMP  CONDITIONS  BY  TWO  APPROACHES 
1.  Eulerian  Coordinates 

Let  i^CilJand  X,('t]be  Eulerian  coordinates  of  the  front  and  back  of 
moving  element  of  mass  which'  is  our  free  body.  Here  is 
independent  of  time.  For  the  observer  fixed  outside  the  flow, 


(/ID 


conservation  of  mass  is  then 

oL  t 

%(t) 

Apply  Leibnitz  rule  (see  last  page) 

Jj  •  J  *  *.  ef***^'  %  *)  -  o 

Note  that  and  X,  in  A  can  be  any  two  fixed  points 
chosen  in  the  flow.  ^ 

Consider  a  control  surface  moving  at  a  velocity  with 
extremes  (t9and  ^(t).  Here  again  eU  is  independent  of  time 
(Eulerian)  xYO  i  J  V-'' 

c-’-) 


X 

o 


Now  and  are  two  positions  in  the  flow  and  can  be 
chosen  coincident  wit)i  and  X,  .  The  .Zp.  in  eq.  A3  is  the 
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same  as  that  in  eq.  A2.  and  represents  the  rate  seen  by  the 
observer  outside  the  flow. 


Combining  eq.  A2  and  A3  by  eliminating  the  integral  with 


gives 


p", 

J 


xYi) 

<*/  f 

‘u  •'VCil 


j(u-U)fJ  -o 


Note  that  and  Also  x, 

are  fixed  points  coinciding  with  yj'  x  "  • 


From  eq.  A4 


_  Now  if  as  the  observer*  s  control  surface  advances,  the 
situation  within  is  independent  of  time  which  is  true  for  a 


detonation  or  stable  deflagration,  then 

This  is  the  condition 


which  is 


the  jump  condition  for  conservation  of  mass,  (eq  18  of  report) 
Now  if  we  rewrite  equation  (Al)  as  (Ala) 

(/{U) 

at 

Then  all  follows  as  previously  leading  to  the  corresponding 


jump  conditions 


and  if,  we  had  used  as  our  physical  statement 

^  *  LCV’  -o 


(^tb) 


then  the  jump  condition  is 


£(^-u)r  f  c  ]  '  =0 


('/IS'iJ 


For  conservation  of  momentum,  the  differential  equation 


corresponding  to  eq  Alb  is 


i/  r  ^ 


The  jump  condition  for  momentum  is  therefore 

2^)  -f  ~0  r  -4  ^ 

Then  using  2^  constant  and  conservation  of  mass  there  follows  the 


jump  condition  for  momentum 


r  - 

[ov  -t  =  0 


(eq  19  of  report) 


For  conservation  of  energy  the  differential  equation 
corresponding  to  eq  Alb  is 

i  a  >',(■<  ')  y, 

The  jump  condition  for  energy  is  therefore 

[(u-U)  ^ 

Using  ^  constant  and  the  jump  condition  for  mass  and  momentum 


this  reduces  to 
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J.  is  the  heat  added  and  is  the  negative  of  the  enthalpy 
constant  h^,  related  to  formation.  In  the  report  if 
represented  h-t  •  jump  condition  equation  19  was 


II.  Lagranaian  Coordinates 

In  these  coordinates,  tags  are  associated  with  certain 
particle  planes.  Conservation  mass  holds  for  any  integration 
limits  in  these  coordinates,  g/y  is  a  function  of  time  s  <r* 


®  f  2.  r"'F4 

^  ’'/o  written  Pf 

•  I'f  "■  •  (■  S  f? '■*)  =£<'■#<?)*! 
=  Aftf-fjf-'-'-f'-f') 

but  ■*  "V' f  ^  -  O  by  cons,  of  mass. 

At  a  particular  time  ^  ^  fj  w  ^  y" 

Using  eq  A3  for  this  integral  gives  * 

il.'"-  ■  £'Ir"  -  *  ^'"•5 


by  cons,  of  mass. 


This  corresponds  to  eq.  A4  and  the  jump  condition  follows  in 


similar  fashion  fo 


added. 
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LEIBNITZ  RULE 


Let  R  be  the  rectangle  defined  by  i  ^  t  i  X' <  ^ 

Let  be  bounded  in  R  and  integrable  with  respect  to  x 

for  each  value  of  t  ,  and  let  ^  exist  in  R  and  also  be 
integrable  in  X  for  each  value  of  ^  .  Let  ^  be 
continuous  for  uniformly  as  to  jf  in  Let  ^/-t) 

be  differentiable  and  y)  continuous  in  x  alone  for  each 

value  of  t  considered  and  for  ^  y  f 


€l(t) 


^here  a  and  b  are  held 
constant,  i.e.,  treated 
as  regular  deMnite  inte¬ 
gral  between  two  fixed 
points. 


dMi) 

clt 


Cl 
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Appendix  B 


SOLUTION  OF  THE  EQUATIONS  OF  THE  ONE  PARAMETER  HYDRODYNAMIC 

COMBUSTION  THEORY 


The  equations  originally  presented  in  the  OPHCT  section  of 
this  report  will  be  repeated  here  for  convenience  when  required  with 
the  original  numbers  retained.  Equations  developed  in  this  appendix 
have  a  letter  B  in  their  designation. 


STEP  1  Solve  eq  32,  33,  34,  35  for  ^ 

in  terms  of  i,  x  and  the  functions  of  temperature  a,  b,  c,  d. 
O 


(3i) 

V  V  f  * 
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'i  “i  -i  -/ 

h  *0  *  9  *0  -  ^k  (^^) 


«0  ■  "o  '  S^  1>  '  ?  A 


Cbj) 


-X 


O  -  /a  0 

j 


We  have  used  ^  ^o)  as  a  change  in  variable  which 

aliminates  fractional  powers  and  gives  a  set  linear  in  x  '■'^  yf)"^ 

Note  that  Vo-C.-j  *'’**  ^  "®‘  appear.  These  will  be 

developed  later. 


STEP  2  Eliminate  x 


combining  eq  26  and  29 


di  Gto  (iiji  (<■) 


X 

duo 


*  \<>  ^  ’b  * 


ft 


v^3 


<2^  /Wo  O  O^ 


H 
-  <r  >> 


Ti?) 


~^-.0  (Bs) 


Eliminate  x^^^  by  combining  eq.  24  and  30. 


1  =  X 


(U. 


f  •*■  X. 


''^o  -  -0/^.  *  'b  ^ 


f<f*i  • 


/V(y<'v 


^  ^  *  *0  *  ’h^  * 


nifOO. 


t  =  / 
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Eliminate  ^  by  combining  equations  B5  and  B6 
Use  f  ®  ^  ^  '’/tdOyi  tor  convenience  in  writing. 

STEP  3-  Substitute  the  results  of  Step  1  in  equation  B7  above. 

This  leads  to  the  following  equation  in  ^  and. 

^  I  C^*0l‘c"g^ 

■t  (it) 

Dividing  thru  by  'f  and  writing  oi  ■k  ^  ir (to-*. 

and  using  S*  *”  *  C"  /**  and  solving  foi 


JL.  - 


'  r/>* 


we  get 


and  solving  for 


*  Since  eq  B9  expresses  in  terms  of  'ST  1  it  can  be  used  in 
eq  21  to  eliminate  .  First,  we  modify  equation  21. 

--  i(K-h)(t  *t) 

STEP  4  -  Use  eq  22  and  23  in  eq  21  and  multiply  thru  by 

^  -i — 

2'x:^ 

^4  ^40,  -»■  J  / ^ 

^  /»».  at  A)  ^  ‘  io,:  ^ 

^  Jt,  O’*  -  r  X,  w,; Z' Z,  V ) 


Z’2.2,) 
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where 


Now  multiply  both  sides  of  BIO  by  since  depends  on  ^ 

only  for  ^  /fy  <“/  B  l  •/. 

.<7; .  ^  (Bi/) 

'  y‘  o/^o-  f 

Since  eq.  B9  gives  f>^  in  terms  of  ^ we  lack  only  expressions 

for  >,  and  x  to  be  able  to  write  an  equation  in  CT  . 

/(a  aijpj  ^  3 

STEP  S-Use  eq  24  to  provide  an  equation  for  ,  multiply  by  ^ 

and  insert  eq 

/>  X  s 

To  obtain  h  v  ,  solve  for  w  - ^  by  combining  eq  28  and  29. 
Use  eq  B6  to  replace  Finally  multiply  by  ^  and  use  eq  B9 


Bt 

.  Bl,  B2,  BS^and  B9  to  obtain  the  equation  for 

V- V”-  wV’' 


fL^/rty 

r  (b/3.) 


This  results  in 


V ''  *  ‘'I  *  V  ’b 

^  p  *’(1/ )/ ^  *^nno>t)  3 

STEP  6- To  obtain  the  polynomial  in  ^ 1  we  substitute  B12,  B14, 
Bl,  B2,  B3,  B4  and  B9  in  Bll,  and  collect  coefficients  of  each 
power  of  ^  .  The  equation  is 


(B>i) 


f?.  B.  'Q  r'*'  CBfs) 
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where 


»e.  A 


II 

f  r/- 

£  =  -^  ( 1  ■*'  ‘^(k)  fi 

A  -  - 

m 

A.''  ♦ 

c. 

^  Y-  Ar  \  rt .  i 
e.  ‘ 

II 

»»»a^  -i 
- c  -#■ 

e, 

v  ' 

n 

VJ^ 

^ac'  -tfr*  PYoia  ('‘“- 

e- 

It 

O  f  1 

4.  = 

/h*  -^ 

■ 

0 

0  ^ 

4 

X?  ore  obtained  from 

4 

K  c. 

respectively  by 

2  -, 

replacing  each  maa./denaity  quotient  by  the  corresponding  !  i 

Thus 
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The  roots  of  eq.  B15  are  more  easily  found  if  it  is  put  in  the  standard 
polynomial  form.  To  achieve  this,  multiply  through  hy  ^  ,  evaluate  the 
product  on  the  left  side,  transpose  and  collect  coefficlente  of  the  Mae 
powers  of  ^  •  This  results  in 

e..  C?"-  -25?^  C  -A?7  ^ 

which  leads  to  the  following 

fi~4)  g'-  ■> 

+  <;,4.)  c'^  ^ (Kg  ■^C.\*K7i'% 

■•■(K  \  B.  -  B-  A  -  6*)  c"+  A  +  C.  <; 

+  (2.,  J>  f  a,?.  +  -Kf>.-C^)g'^ 

+  r,^,)  5-' 

+  ('d.f.  zi,^.  - 

+  (j>^  f.  ^  c-p )  =  o  J 
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Equation  Bl6  is  the  required  polynomlnal  in  ^  vhlcb  conblneB  all  10 
equations  and  avoids  the  need  for  iterative  solutions  of  the  set.  Once  ^ 
is  found,  the  various  equations  are  used  to  obtain  for  a  particular 
the  corresponding  values  of  fa  e-  Xhen  ^2  follows  from  eq. 

22.  The  velocities  ^  are  obtained  from  the  Juaqp  conditions 

for  aass  and  acBentum. 
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